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Abstract 

The  feedout  process  transfers  mass  perturbations  from  the  rear  to  the  front  surface  of  a 
driven  target,  producing  the  seed  for  the  Rayleigh-Taylor  (RT)  instability  growth.  The  feedout 
mechanism  is  investigated  analytically  and  numerically  for  the  case  of  perturbation  wavelength 
comparable  to  or  less  than  the  shock-compressed  target  thickness.  The  lateral  mass  flow  in  the 
target  leads  to  oscillations  of  the  initial  mass  non-uniformity  before  the  reflected  rippled 
rarefaction  wave  breaks  out,  which  may  result  in  RT  bubbles  produced  at  locations  where  the 
areal  mass  was  initially  higher.  This  process  is  determined  by  the  evolution  of  hydrodynamic 
perturbations  in  the  rippled  rarefaction  wave,  which  is  not  the  same  as  the  Richtmyer-Meshkov 
(RM)  interfacial  instability.  An  exact  analytical  formula  is  derived  for  the  time-dependent  mass 
variation  in  a  rippled  rarefaction  wave,  and  explicit  estimates  are  given  for  the  time  of  first  phase 
reversal  and  frequency  of  the  oscillations.  The  limiting  transition  from  the  case  of  RM 
perturbation  growth  at  large  density  difference  (low  ambient  density  behind  the  rear  surface)  to 
the  case  of  feedout  (zero  density)  is  studied,  and  it  is  shown  that  the  latter  limit  is  approached 
only  if  the  ambient  density  is  extremely  low,  less  than  1/1000  of  the  pre-shock  target  density. 
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I.  INTRODUCTION 


One  of  the  sources  that  seed  the  initial  mass  perturbations  for  the  Rayleigh-Taylor 
(RT)  perturbation  growth  in  laser  fusion  targets  is  the  roughness  of  the  inner  surface  of 
the  target.  The  process  that  generates  the  RT  mass  perturbation  seeds  from  the  initial 
inner  surface  ripples  is  called  feedout,  see  Refs.  1  to  3  and  references  therein.  It  should  be 
noted  that  the  first  discussion  and  numerical  simulation  of  the  inner  target  surface 
roughness  seeding  the  RT  instability  at  its  outer  surface  were  presented  in  Ref.  4.  In  this 
paper,  we  neglect  all  other  sources  of  non-uniformity  (laser  imprint,  front  surface 
roughness)  and  assume  a  constant  pressure  drive  that  may  correspond  either  to  a  low- 
energy  foot  of  the  laser  pulse  or  to  a  constant  hohlraum  temperature.  Then  for  the  single¬ 
mode,  long-wavelength  case1'3  defined  in  planar  geometry  by 

X»2nLs,  (1) 

where  A  is  the  perturbation  wavelength,  and  Ls  is  the  post-shock  target  thickness,  the 

feedout  could  be  visualized  as  shown  in  Fig.  1 .  A  planar  shock  wave  S  approaches  the 
rear  surface  ripples,  first  breaking  out  at  the  valleys,  and  then  at  the  peaks  [Fig.  l(a,  b)]. 

A  rippled  rarefaction  wave  is  reflected  first  from  the  valleys,  then  from  the  peaks.  In  the 
reference  frame  of  the  shocked  mass,  its  leading  edge  L  propagates  at  constant  velocity 
«, ,  the  speed  of  sound  in  the  unperturbed  shocked  material.  Therefore,  the  leading  edge 
of  the  reflected  rarefaction  wave  is  a  sine  wave  of  the  same  phase  as  the  initial  rear 
surface  ripple,  and  time-independent  amplitude  [Fig.  1(c)].  The  shocked  material  expands 
into  vacuum,  starting  from  the  valleys,  at  a  constant  velocity  ue ,  which  for  ideal  gas 
equals  2a,  /(y,  - 1) ,  so  that  after  the  sine-shaped  expansion  front  E  is  fully  formed,  its 
ripple  amplitude  is  also  time-independent.  The  expansion  front  propagates  faster  than  the 
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incident  strong  shock  wave  in  the  target,  if  y1  <  2  +  y/5  =  4.236,  in  which  case  the  phase 
of  the  expansion  front  is  inverted  with  respect  to  the  initial  rear  surface  ripple  [Fig.  1(c)], 
Breakout  of  the  rippled  rarefaction  wave  occurs  when  its  valleys  (originated  from  the 
valleys  at  the  rear  surface  of  the  target)  reach  the  front  surface  [Fig.  1(c)],  starting  its 
acceleration.  In  other  words,  the  thinner  parts  of  the  target  start  accelerating  earlier,  and 
experience  a  higher  acceleration  under  the  same  driving  pressure  because  of  the  lower 
areal  mass.  During  the  RT  growth  phase  that  follows,  these  thinner  parts  evolve  into 
bubbles,  propagating  ahead  and  dumping  more  of  their  mass  into  the  spikes  that  trail 
behind.  This  physical  picture  has  been  recently  studied  theoretically1  and 
experimentally  . 

In  the  present  paper  we  study  the  case  when  the  long-wavelength  condition  (1) 
does  not  hold.  Then  the  situation  changes  in  one  important  way:  the  lateral  mass 
redistribution  in  the  rippled  rarefaction  wave  can  no  longer  be  neglected.  The  lateral  mass 
flow  is  driven  as  shown  in  Fig.  1(b).  After  the  shock  breaks  out,  the  expansion  starts  from 
the  valleys,  decreasing  the  pressure  there,  whereas  near  the  peak  locations  the  pressure  is 
maintained  at  the  constant  post-shock  value.  The  resulting  lateral  pressure  gradient  starts 
driving  the  mass  from  the  peaks  to  the  valleys,  decreasing  the  pressure  at  the  peaks  and 
increasing  it  at  the  valleys.  The  mass  flow  continues  when  the  pressure  gradient  vanishes, 
thus  overshooting  the  equilibrium  situation,  and  building  up  a  reversed  pressure  gradient. 
These  sonic  oscillations  of  areal  mass  in  a  rippled  rarefaction  wave,  which  continue  until 
the  rarefaction  wave  breaks  out  at  the  front  surface  of  the  target,  were  first  described  in 

Ref.  5  (Fig.  14).  In  the  context  of  the  feedout  problem,  they  were  observed  in  the 
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hohlraum-driven  experiments  on  Nova  and  reproduced  in  the  simulations. 
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Feedout  as  described  above  is  a  hydrodynamic  process  that  generates  the  actual 
seed  for  the  RT  instability  starting  from  the  ripples  at  the  rear  surface  of  the  target.  It 
should  be  emphasized  that  the  feedout  mechanism  per  se  does  not  involve  any  interfacial 
hydrodynamic  instability.  There  is  no  perturbation  growth  neither  at  the  leading  edge  of 
the  rippled  rarefaction  wave  nor  at  the  rippled  expansion  front.  The  density  and  pressure 
perturbations  within  the  rippled  rarefaction  wave  evolve  as  decaying  sonic  waves.  Note, 
however,  that  the  stable  feedout  situation  could  be  regarded  as  a  limiting  case  of  a  RM 
unstable  flow,  if  we  assume  some  low  ambient  density  p2  behind  the  rear  surface  of  the 
target,  and  then  let  p2  tend  to  zero.  For  any  small  but  finite  p2  the  rippled  rarefaction 
wave  is  unstable,5' 6  and  the  RM  instability  growth,  asymptotically  linear  in  time,  must 
take  place  at  the  contact  interface  between  the  decompressed  target  material  and  the 
shock-compressed  ambient  gas.  Manifestation  of  the  RM  instability  has  been  detected  in 
the  feedout  simulations  of  Ref.  2.  Do  we  have  a  contradiction  here?  To  answer  this 
question,  it  is  important  to  study  in  detail  the  limiting  transition  p2  — >  0  from  the  RM 
unstable  to  a  stable  feedout  situation. 

This  paper  is  structured  as  follows.  In  Section  II,  we  describe  analytically  the 
evolution  of  mass  perturbations  in  a  rippled  rarefaction  wave  produced  when  a  planar 
shock  wave  hits  a  rippled  rear  surface  of  the  target,  compare  the  results  to  the  RM 
unstable  situation  when  the  half-space  beyond  the  rear  surface  of  the  target  is  filled  with  a 
low-density  gas,  and  study  the  limiting  transition  between  the  two  cases.  In  Section  III, 
we  present  simulation  results  for  both  cases  and  discuss  the  opportunities  for  their 
experimental  verification,  and  Section  IV  concludes  with  a  discussion. 
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II.  THEORY 


A.  Formation  of  a  rippled  rarefaction  wave 

Consider  a  planar  layer  of  ideal  gas  whose  adiabatic  exponent,  initial  density,  and 
thickness  are  yl ,  p01 ,  and  L0 ,  respectively.  The  half-space  behind  the  layer  x>  L0 

contains  a  low-density  gas  with  density  p2  and  adiabatic  exponent  y2 .  We  are  interested 
in  the  case  of  low  ambient  density  p2  «  p01  and  in  the  limiting  transition  to  the  case  of 


feedout,  p2  =  0  .  To  exclude  expansion  of  the  layer  into  vacuum  from  behind  prior  to  the 
shock  arrival  in  the  latter  case,  we  neglect  the  initial  pressure  in  both  gases  compared  to 
the  shock  pressure.  Then  a  constant  pressure  px  instantly  applied  to  the  front  layer 


surface  produces  a  strong  shock  wave  whose  velocity  is 


(Yi  +  DPi 

^Poi 


(2) 


Shock  compression  equals  (7,  +  1  )/(7,  - 1) .  Density,  x-velocity  and  speed  of  sound  in  the 
shocked  gas  are 


- rPoi>  vo  - r^o»  a\ 

Yi  “I  7i  +  1 


[27,(7, -I)]' 


respectively.  By  the  end  of  the  shock  transit  time  interval,  ts  =  L0/  D0,  the  layer  is 
compressed  to  its  post-shock  thickness,  Ls  =  (y1  - 1  )L0  l((yx  + 1) . 

At  the  instant  of  shock  breakout,  a  centered  rarefaction  wave  starts  to  propagate 
into  the  uniform  shocked  gas.  Let  us  pass  to  the  reference  frame  moving  with  the  shocked 
gas  and  make  t  =  0  and  x  =  0  correspond  to  the  time  and  place  of  the  shock  breakout, 
respectively.  The  gas  flow  in  the  centered  rarefaction  wave  is  self-similar: 
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(4) 


A(x,t) 


a(x,t) 

a. 


2 

7i  +  l 


7 i-l  * 
Yi+l  axt  ’ 


p(x,t)  =  P\A{x,t)2/i7l'l) ,  p(x,t)  =  plA(x,t)2nKn'1)  ,  (5) 

where  a(x,t) ,  p(x,t  )  and  p(x,t  )  are  local  values  of  the  sound  velocity,  density,  and 
pressure,  respectively.  The  strength  of  the  rarefaction  wave  is  characterized  by  the  Mach 
number 


4T,  =  — ,  (6) 

ax 

where  ax  is  the  speed  of  sound  behind  the  rarefaction  wave.  The  normalized  speed  of 
sound,  A(x,  t) ,  in  a  rarefaction  wave  varies  between  1  at  its  leading  edge  and  Mx  at  its 
trailing  edge.  The  values  of  density  and  pressure  at  the  trailing  edge  are  found  by 
substituting  A- Mx  into  (5): 

Pi-  p\=  pxMxr,Kr'~l).  (7) 

The  Mach  number  Mx  varies  between  1  and  0,  the  former  value  corresponding  to  a  weak 
rarefaction  (sound  wave),  and  the  latter  —  to  a  strong  rarefaction  (expansion  into  vacuum). 
In  our  case  of  strong  shock-interface  interaction,  its  value  is  found  from 

_P 2_  =  Oh+1  M  (8) 

Poi  (y2+l)[  +  s1/2(l-M1)J’ 

where  s  =  2y,  /( y,  - 1) ,  see  Refs.  5,  7.  Inspection  of  Eq.  (8)  reveals  that  expansion  into 
vacuum  (Mx  =  0 )  corresponds  to  the  feedout  case  p2  =  0.  However,  a  small  value  of  the 
ratio  p2  /  p01  does  not  necessarily  imply  that  Mx  is  also  small  because  the  power  s  could 
be  large.  Taking,  for  example,  7\  ~  7 2  =  5  /  3  ,  so  that  s  =  5 ,  we  find  that  density  ratios 
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p 2  /  p01  =  10  10  2. 10  and  1 0  4  would  produce  Mach  numbers  of  the  rarefaction  wave 

Mj  =  0.75,  0.53,  0.36  and  0.24,  respectively.  Note  that  the  self-similar  profiles  (4),  (5) 
do  not  depend  on  the  rarefaction  Mach  number  M, . 

The  leading  edge  of  the  rarefaction  wave  propagates  in  the  negative  x  direction 
into  the  unperturbed  resting  uniform  fluid  at  local  sound  velocity,  cp :  x,(t)  =  -axt  .  The 


rarefaction  transit  time  equals 


(9) 


where  470  is  the  Mach  number  characteristic  of  the  shocked  fluid  flow  with  respect  to  the 
shock  front.  For  an  arbitrary  equation  of  state,  it  must  be  M0  <  1 ,  so  that  the  rarefaction 
transit  time  is  always  less  than  the  shock  transit  time. 

The  trailing  edge,  which  also  propagates  at  local  sound  velocity  a]  with  respect  to 
the  fluid  particles,  separates  the  rarefaction  wave  from  the  uniform  downstream  flow, 
where  density  and  pressure  have  constant  values,  p{  and  p* ,  and  the  fluid  velocity 
equals 

=  (1_Ml)’  (10) 

y  i  i 

so  that  xt(t)  =  (ut  -Mxax)t.  For  expansion  into  vacuum,  0 ,  we  have:  p*  =  0  , 
p\-  0  ,  there  is  no  gas  behind  the  trailing  edge,  which  in  this  case  is  called  the  expansion 
front.  Its  velocity  equals  ue  =  2«,  !(yx  - 1)  and  coincides  with  the  fluid  particle  velocity 
(10). 
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Now  introduce  a  small  single-mode  ripple  at  the  rear  surface  of  the  layer. 
Denoting  the  initial  ripple  amplitude  by  Sx0 ,  we  find  that  immediately  after  the  shock 
breakout  the  ripple  amplitudes  at  the  leading  and  trailing  edges  of  the  rarefaction  wave 
are7 


5x]  =  — - 1  +  [2 y,  (7!  - 1)]  2  >*0  , 

7i  +1 


sx  ;  =  — 7.-1+ 

7i  +  l  7,-1 


[(7i+1)Mi_2]Ux0. 


The  leading  edge  ripple  perturbation  amplitude  is  time-independent,  8x,(t)  =  Sx*  .  The 
trailing  edge  perturbation  amplitude  grows  linearly  with  time5’ 6  for  any  nonzero  Mx .  In 


the  case  of  expansion  into  vacuum  (feedout),  it  is  also  time-independent:  5xt(t)  =  Sx* . 

B.  Areal  mass  perturbations  in  a  rippled  rarefaction  wave 

A  small-amplitude  stability  analysis  of  a  rippled  rarefaction  wave  has  been  done 
analytically,  first  in  Ref.  8,  and  later,  in  a  different  way,  in  Ref.  5,  and  numerically,  in 
Ref.  6.  The  results  presented  below  are  based  on  the  analysis  of  Ref.  5,  where  the  details 
could  be  found.  Since  the  density  is  continuous  at  the  edges  of  the  rippled  rarefaction 
wave,  the  contribution  of  the  rippled  rarefaction  wave  to  the  total  areal  mass  perturbation 
(other  contributions  are  discussed  in  Section  IIIC)  is  given  by 


x,(t)  x,(t) 

Smr—  1 8pdx-pl  | 

x,(t)  x,(t) 


Sp(x,t ) 

p(x,t ) 


dx 


=  pxSx*T 


7,+l 

7,(7,  -1) 


\A1/{7l-')xV(A,x)dA. 

Mj 


(13) 


Here, 


(14) 


T  =  kaxt 


is  normalized  time  ( k  =  2 n  I X  is  the  perturbation  wave  vector),  Sp  is  the  local  density 
perturbation,  is  the  relative  pressure  perturbation  in  the  rarefaction  wave 

normalized  with  respect  to  the  constant  value  k&c* ,  and  A  is  related  to  x,  t  by  (4).  An 
explicit  expression  for  x¥(A,t)  is  available  for  A  close  to  1: 


r  \1/2 

7i+l 

ri-1 


(1  -  A)1/2t 


(15) 


where  J0(z)  is  the  Bessel  function.  In  a  general  case, 


1  ^ 

»p(^t)  =  -£  a 


./• 


Y/i+O 

T2  " 

4 

(16) 


where  calculation  of  the  functions  *Pj(A),  j  -  0, 1, 2,...,  is  described  in  Ref.  5.  In 
particular, 

2y,  (1-13  +  BA) 

VV  •  <17) 

where,  in  our  strong-shock  approximation, 

2y1(y1  +  i) 


(8  = 


2y1(y1-l)  +  (2y1)1,-!(71-l)^ 


(18) 


In  the  limit  A  — >  1  this  is  consistent  with  (15),  as  it  should  be.  With  the  help  of  ( 1 6)-(  1 8), 
we  find  the  initial  value  of  areal  mass  perturbation  characteristic  of  the  rippled  rarefaction 


wave: 


5m r  (°)  _  1  _  7i  ~  1  0  +  M2/(r,-i) 


PAi 


Yi  +1 


(8-1- 


V 


2 /3M, 
7i  + 1 


.  (1  -  M, )  ,, 

2V - 44,  M,  — >  1, 

Yi  -1 


1 


Y\  ~  1 
Yi  + 1 


(19) 


/3,  M,  -»0. 
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A  weak  reflected  rarefaction  wave  could  be  produced  in  our  case  of  a  strong  shock- 
interface  interaction  if  the  density  ratio  p2  /  p01  were  a  little  less  than  the  total 

transmission  value  (Y,  + 1)/(/2  +  0  found  by  substituting  Ml  =  1  into  (8). 

(Alternatively,  the  density  ratio  could  be  much  smaller  than  this  value,  but  then  the 
incident  shock  wave  should  be  weak.)  For  instance,  assuming  Yi  ~  Yi  ~  5/3  and 
Pi !  Poi  =  0-95  ,  we  find  from  (8):  Ml  =  0.9946.  Figure  2  shows  the  time  history  of  areal 

mass  perturbation  calculated  analytically  for  this  case.  The  solid  and  dotted  lines 
correspond  to  the  exact  solution  (20)  and  to  the  approximation  (21),  respectively,  in  good 
agreement  with  each  other.  The  oscillations  are  standing  sonic  waves  excited  in  the 
rarefaction  wave  by  the  lateral  pressure  gradient,  which  drives  the  gas  from  the  high  to 
the  low  pressure/density  areas.  The  lateral  mass  flow  continues  when  the  pressure 
gradient  vanishes,  thus  overshooting  the  equilibrium  situation,  but  not  by  much:  as  seen 
in  Fig.  2,  the  peak  value  of  mass  perturbation  amplitude  with  inverted  phase  is  only  13% 
of  its  initial  value.  At  late  time,  the  oscillation  amplitude  decays  as  t^3/2>  see  Eq.  (21). 
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For  the  case  of  feedout  (strong  rarefaction  wave,  Mx  =  0 ),  the  integrals  in  Eq. 


(20)  are  evaluated  explicitly,  without  calculating  the  functions  Tf(/I) .  This  is  done  by 

multiplying  Eqs.  (63),  (64)  of  Ref.  4,  which  determine  these  functions,  by  A1,(y'  '' , 
integrating  them  over  A  from  0  to  1,  and  using  the  boundary  conditions  (66).  Then,  after 


some  algebra,  we  obtain: 


r  20v-l)  3 

dm  (t)  =  p0I&0 •  2(3-/.)(y.+i)£H)i r liilY] 1  L  3-y. _ 2 

’  01  °  y1+i-(y1-i)i8feTM.[3-r1j4j  r  y1+!  ,4  l 

3 -7/  2 


- l)(2y  -  >3)  +  y,  + 1] 

[(7:  - 1)7  + 1]{4[(7,  - 1)7  +  IT  - (7, -l)2}’ 


where  r(z)  is  the  Euler  gamma  function.  Note  that  this  expression  is  exact,  not 
approximate,  like  (21),  and  its  right-hand  side  is  finite  at  7,  =3.  Equation  (19)  was  used 
to  normalize  the  areal  mass  perturbation  amplitude  with  respect  to  its  initial  value,  which, 
of  course,  is  the  same  as  the  rear  surface  mass  non-uniformity  of  the  target  prior  to  the 
shock  breakout:  Smr  (0)  =  pm&c0 . 


The  infinite  series  (22)  converges  at  all  times  and  represents  an  entire  analytic 
function  of  complex  time  T.  At  early  time,  this  function  could  be  approximated  by  the 
first  two  terms  of  the  series: 


Smr(T)  =  p0l8x0 


at  t  <  fj , 


(23) 


where  the  coefficient 

T  =  / 27i(37i  -  l)l7i+l-(7i  -I)/?]]/'2 
1  (7i+l)[37i-l-(7i-l)0]  J 


(24) 
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could  be  regarded  as  the  first-order  estimate  of  the  time  interval  between  the  shock 
breakout  and  the  first  phase  reversal  of  areal  mass  perturbation. 

At  |  t  |— »  oo  ,  the  function  generated  by  the  series  (22)  has  an  essential  singularity, 
like  exponential  or  Bessel  functions.  The  main  exponential  term  describing  its  divergence 
near  this  singularity,  at  t  -4  ±z°°  ,  is  immediately  found  by  simplifying  the  expression  for 
the  series  coefficients  at  large  j:  Smr  =  const  •  sinh(Qr)/(Or) ,  where  the  values  of  const 
and  Q.  are  given  below.  On  the  real  axis,  the  corresponding  terms  become  oscillatory, 
describing  the  most  rapidly  varying  part  of  the  solution: 


5mr(T)  =  p0lSx0 


4yfl(y  j+1)  sin(Of) 

-l)[7i +l-(7i -0/5]  Gt 


+  (...)  , 


where  ( )  stands  for  the  terms  slowly  varying  in  comparison  with  sin(f>r)/(f>T) , 


'2(7. -1)" 

7i  +1 


(26) 


is  a  function  of  7.  that  decreases  from  1  at  — >  1  to  1/2  at  7i  — >  00  .  In  particular, 

Q(7.  =  5/3)=  1/V2  =  0.707  and  Q(y,  =  3)=  1/ V^=  0.607  ■ 

Figure  3  demonstrates  oscillations  of  areal  mass  calculated  using  (22)  for  strong 
rarefaction  waves  ( M ,  =  0 ,  expansion  into  vacuum),  which  emerge  in  the  feedout 
process.  The  first  term  in  the  right-hand  side  of  Eq.  (25)  turns  out  to  be  a  good 
approximation  to  the  late-time  behavior  of  5mr  for  I  <  7,  <2.  in  Fig.  3(a),  we  compare 
the  exact  solution  (22)  for  7i  =5/3  (solid  line)  to  the  asymptotic  term  (25)  (dotted  line), 
demonstrating  a  good  agreement  after  f  ~  15.  Figure  3(b)  compares  the  exact  results  for 
different  values  of  7i  between  1  and  2.  Overshooting  of  the  initial  amplitude  is  found  to 
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be  large  for  all  cases:  by  a  factor  ranging  from  2.5  to  3.3  for  y,  varied  between  7/5  and 
9/5.  One  or  two  of  the  following  peaks  are  even  higher.  Figure  3(c)  demonstrates  a 
different  behavior  found  for  higher  y  s.  The  amplitude  is  shifted  to  negative  values  and 
exhibits  a  slow  (not  faster  than  linear)  growth,  which  becomes  more  pronounced  as  y  is 
increased.  This  could  indicate  some  weak  instability  of  a  rippled  rarefaction  wave  that 
takes  place  even  at  M{  =  0 ,  is  visible  mostly  at  high  y  s,  and  has  not  been  noticed 
before.5’ 6  Note  that  the  oscillating  contributions  to  the  exact  solutions  shown  in  Fig.  3(c) 
are  well  approximated  by  Eq.  (25). 

The  time  of  first  phase  reversal  is  in  good  agreement  with  the  estimate  (24),  the 
error  in  all  cases  being  4%  to  5%.  In  agreement  with  (26),  frequency  of  oscillations 
slightly  decreases  with  increased  Y,  ■  At  late  time,  the  oscillation  amplitude  decays  as 

t  \  as  predicted  by  Eq.  (25).  The  time  interval  between  the  first  and  the  second  phase 
reversals  scales  as  1  /  £2  : 

At,  =^e(Yl),  (27) 

where  0(7,)  is  a  dimensionless  correction  factor  of  order  unity.  This  factor  tends  to  unity 
in  the  limit  y,  — >  1 .  Its  values  for  7i  =  7/5  ,  5/3,  9/5,  2,  8/3,  and  2.82  are  1.04,  1.19,  1.26, 
1.36,  1.77,  and  2.03,  respectively.  As  seen  in  Fig.  3(c),  for  7,  >  2.82  ,  no  second  phase 
reversal  is  predicted. 

Figure  4  shows  how  the  mass  perturbation  behavior  changes  as  its  Mach  number 
is  gradually  increased  from  0  to  0.8  (7i  =  Yi  =  5  /3  ).  As  found  from  (8),  the  values  of 
M i  =  0.2,  0.4,  0.6,  and  0.8  correspond  to  the  ratios  p2  /  p0l  =4.17E10"5,  1.97E10"3,  2.27E10"2, 
and  0.16,  respectively.  To  facilitate  the  comparison,  the  mass  perturbation  amplitudes  are 
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also  normalized  here  to  the  respective  values  of  Smr(0)  given  by  (19).  We  see  that  the 
amplitude  of  the  first  inverted  peak  is  very  sensitive  to  the  density  ratio.  E.  g.,  addition  of 
a  very  small  ambient  density  p2  /  p01  =  1.9/E10  3  behind  the  rear  surface  results  in  a 
significant  (27%)  reduction  of  the  relative  amplitude  of  the  negative  peak  compared  to 
the  case  of  feedout  (expansion  into  vacuum). 

This  theory  of  feedout  applies  until  the  rarefaction  wave  breaks  out  at  the  front 
surface.  From  (9)  and  (14),  we  find  that  this  instant  corresponds  to  z  =  zr  defined  as 


zr  =  kaxtr  =  kLs  = 


A  ’ 


(28) 


cf.  Eq.  (1).  If  zr  >  T) ,  then  our  feedout  experiment  or  simulation  must  show  a  phase 
reversal  of  mass  perturbation  in  a  rippled  rarefaction  wave.  If  zr  >  r,  +  n0  i  20 ,  then  we 
should  also  observe  the  first  inverted  peak.  For  the  experiments  of  Ref.  3  (plastic,  L0  =25 
pm,  A  =  100  pm),  assuming  a  (2.5  to  3)-fold  shock  compression  of  plastic,  we  find: 
zr  =  0.52  to  0.63  .  The  estimate  (24)  demonstrates  that  this  is  definitely  insufficient  for 
phase  inversion:  in  the  range  of  Y\  from  1  to  2.5  the  value  of  t,  varies  between  1.25  and 
1.64.  As  found  from  (23),  the  above  value  of  f,.  corresponds  only  to  a  ~1 5  to  20% 
reduction  in  the  mass  variation  amplitude  by  the  time  of  the  rarefaction  wave  breakout  in 
comparison  with  the  long-wavelength  limit  zr  — >  0 .  Thus  the  long-wavelength 
assumption  (1)  is  valid  in  the  case  of  this  experiment,  and  the  theory  of  Ref.  1  applies. 
This  is  not  the  case  for  the  thick-foil  experiment  of  Ref.  2,  where  the  initial  thickness  of 
an  aluminum  slab  was  L0  =  86  pm,  and  the  rear  side  ripple  wavelength  A  =  50  pm.  A  3- 
fold  shock  compression  of  solid  aluminum  would  correspond  to  T  =  3.6  ,  which  exceeds 
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Tj  for  any  reasonable  value  of  yl .  In  particular,  for  y1  =  5/3  we  have: 

Tj  +7i0 / 20  =  4.0  .  As  demonstrated  by  Eq.  (24)  and  our  Fig.  3,  the  value  of  zr  =  3.6  is 
large  enough  for  phase  inversion  and  agrees  with  a  ~2. 5-fold  excess  of  mass  variation 
with  inverted  phase  over  its  initial  value,  which  is  shown  in  Fig.  7  of  Ref.  2. 

C.  RM  instability  at  large  density  difference 

We  have  demonstrated  above  that  the  feedout  process  does  not  involve  any 
interfacial  hydrodynamic  instability.  Nevertheless,  the  vorticity-driven  RM  perturbation 
growth  at  the  rear  surface  has  been  discussed  in  connection  with  feedout  in  Refs.  2,  3.  To 
find  out  if  there  is  any  contradiction  here,  let  us  assume,  as  done  in  most  computational 
studies  that  use  Eulerian  or  Fagrangian-Eulerian  codes,  that  there  is  some  low-density 
ambient  gas  behind  the  rear  surface  of  our  gas  layer,  and  consider  the  limiting  transition 
from  the  low-density  case  p2  /  p01  «  1  to  the  feedout  case  of  p2  -  0.  In  the  former  case, 

the  RM  instability  would  inevitably  develop  for  any  p2  /  p01,  no  matter  how  small, 
whereas  the  latter  case  is  essentially  stable.  Obviously,  we  must  obtain  the  same  results 
for  evolution  of  mass  perturbations  in  the  stable  case  p2  =  0  and  in  the  limiting 
transition  p2  /  p01  — >  0  from  an  essentially  RM  unstable  situation.  The  questions  are: 

How  small  should  the  ratio  p2  /  p01  be  in  order  to  make  all  the  contributions  to  areal 

mass  variation,  other  than  the  rippled  rarefaction  wave,  negligible?  How  does  the  RM 
instability  develop  when  its  contribution  to  the  mass  variation  is  no  longer  significant? 

The  shock-interface  interaction  makes  the  contact  interface  move  in  positive  x- 
direction  at  the  velocity  ut  given  by  (10)  and  produces  a  strong  transmitted  shock  wave 
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propagating  into  the  ambient  gas  at  the  velocity  D2  =  (y2  + 1  )ut  /  2  .  Now  the  contact 
interface  separates  ambient  gas,  shock-compressed  to  the  density 

p2=— \p2,  (29) 

72-1 

and  the  expanded  rear  surface  of  our  layer,  where  the  gas  density  is  given  by  (7).  From 
(7),  (8)  and  (29)  we  find  the  post-interaction  density  ratio  at  the  contact  interface: 

P;_  ov -i)Mt2  _  Ob  - 1) 

Pi  (y2  -1)|  +  51/2(1-A/'1)]2  (r2 -!)(!  +  51/2 )2/r‘ 

in  the  limit  p2  /  p0l  — >  0  .  The  density  contrast  at  the  interface  therefore  remains  high  in 

this  limit,  although  the  densities  at  both  sides  of  the  interface  tend  to  zero.  At  p2  =  0  ,  the 

contact  interface  coincides  with  the  trailing  edge  of  the  rarefaction  wave,  the  expansion 

front,  see  (10). 

The  total  mass  variation  produced  in  the  RM  unstable  flow  could  be  presented  as 
5m(z)  =  Smr  (f)  +  8mse (f)  +  8mc (r)  +  5m  (t).  (3 1 ) 

Here,  8mr(z)  is  the  mass  variation  in  the  reflected  rarefaction  wave  given  by  (20); 

5mse(r)  is  the  contribution  due  to  sonic  and  entropy  perturbations  in  the  uniform  flow 
regions  on  both  sides  of  the  contact  interface, 

Smjz)  =  1 8p(x,z)  dx,  (32) 

[where  5p{x,z  )  is  the  local  density  perturbation,  and  xs  =  D  ,t  is  the  coordinate  of  the 
transmitted  shock  front],  8mc(z)  and  8ms(z)  are  contributions  of  the  perturbations  at  the 
contact  interface  and  transmitted  shock  front,  respectively: 


(72  +1)P2 


nvn-u'  7i 


(7i+l)Po 


(30) 
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8mc (t)  =  (pi  -pl)  Sxc (t),  Sms (t)  =  (p2  -  p2 )  (T), 


(33) 


[where  8xc(z)  and  8xjz)  are  the  displacement  amplitudes  of  the  contact  interface  and 
transmitted  shock  front,  respectively].  We  calculate  the  time-dependent  perturbation 
amplitudes 8p(x,z),  8xc(z),  and  Sxjz),  and  evaluate  8m(z)  analytically  following  the 
procedures  of  Refs.  7,  9,  and  using  Eq.  (20).  The  results  shown  in  Figs.  5(a)  to  (d)  refer  to 
Yi  -  7  2  =  5/3  and p2  /  p01  =  10  1 ,  1 0  2 ,  10  and  1 0  4 .  We  plot  the  total  mass  variation 
(31)  and  the  contributions  to  it  due  to  rippled  rarefaction  wave,  8mr(z)  (which  should 
dominate  in  the  feedout  limit)  and  to  the  rippled  contact  interface,  8mc(z )  (this  one 
dominates  at  late  times  in  the  RM-unstable  situation). 

For  p,  /  p01  =  10  1  [Fig.  5(a)]  the  linearly  growing  interfacial  contribution  8mc(z ) 
to  the  total  mass  variation  is  dominant,  whereas  the  lateral  mass  transfer  in  the  rippled 
rarefaction  wave  adds  some  oscillatory  correction.  For  p2  /  p01  =  10  2  [Fig.  5(b)], 
contribution  of  the  oscillations  is  much  larger,  but  the  linearly  growing  term  Smc(z )  still 
prevails  in  determining  the  evolution.  At  p2  /  p01  =  10  3  [Fig.  5(c)],  the  contribution  of 
the  rippled  rarefaction  wave  characteristic  of  the  feedout  is  dominant,  although  the 
correction  due  to  the  interfacial  term  is  not  so  small  —  about  1 0%  at  the  first  inverted 
peak.  Finally,  at  p2  /  p01  =  |  (F4  [Fig.  5(d)]  the  contribution  of  the  rippled  rarefaction  wave 

prevails.  We  come  to  the  conclusion  that  an  accurate  numerical  modeling  of  the  feedout 
process  requires  the  ambient  density  (which  is  not  there  in  the  feedout  experiments,  and 
is  typically  introduced  in  simulations  to  make  the  codes  work)  to  be  less  than  10"3. 
Otherwise  the  simulation  results  could  contain  a  numerical  artifact  —  that  is,  a  significant 
contribution  from  the  contact  interface  which  is  not  really  there. 
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We  have  seen  that  in  the  limit  p2  /  p0l  — >  0  the  contribution  of  the  contact 
interface  to  the  total  mass  variation  becomes  small.  Its  smallness,  however,  comes  from 
the  small  density  factor  p~  -  p2  in  Eq.  (33).  It  does  not  imply  that  the  interface 
displacement,  8xc( t),  should  also  be  small.  In  fact,  8xc (t)  typically  exceeds  the  initial 
value,  Sx0 ,  by  a  large  factor,  growing  with  time  as  it  should  in  the  RM  unstable  situation. 
This  perturbation  development,  as  we  demonstrate  below,  is  somewhat  unusual  due  to  the 
low  value  of  p2  /  p0l.  Besides,  in  the  low-density  limit  we  can  expect  a  freeze-out10, 
which  means  zero  asymptotic  value  of  the  RM  interfacial  growth  rate.  For  our  case  of 
reflected  rarefaction  wave,  the  latter  could  be  estimated  with  the  aid  of  the  so-called 
prescriptions  due  to  either  Meyer  and  Blewett  (MB)11  or  Vandenboomaerde  et  al. 
(VMG)12: 


lim — 8x(t)  = 


kut  x 
ku,  x 


A T  (&r0  +  8x*c )  (MB ) , 

'  (A5y0  +AS&;  )  (VMG), 


(34) 


where  the  pre-  and  post-interaction  Atwood  numbers  are  defined  as 
A  =  (p2  -  p01)/(p2  +  p01)  and  A *  =  (p*  -  p*)/((p2  +  p*) ,  respectively,  and  Sx*  is  the 
post-interaction  interface  ripple  amplitude,  which  in  our  strong-shock  case  equals  to 

1 


Sx.  = 


7  i+l 


y1-l-2(l-M1j 


{  2 yt  ) 

1/2  j 

y,  -1 

V1  J 

] 

Sxn 


(35) 


[at  M j  =  0 ,  the  contact  interface  coincides  with  the  expansion  front,  so  that  the 
displacement  amplitudes  (35)  and  (12)  are  identical,  as  they  should  be.]  In  the  limit 
p2  /  p01  — >  0  both  A  — »  -1  and  A*  — >  -1 ,  so  both  prescriptions  predict  nearly  the  same 
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asymptotic  growth  rate.  The  freeze-out,  r„  =  0  ,  is  predicted  for  the  case  when  the  shock- 
interface  interaction  inverts  the  phase  of  the  interface  ripple  without  affecting  its 
amplitude,  8x *  =  -8x0 .  The  corresponding  density  ratio  is  found  from  (8),  (34),  (35): 


P2_ 

.  freeze-out 


(r,  +  i)(r2  +  o 


7i(7i  -1) 

2 


1/2  'l  27l/(ri-D 


The  RM  freeze-out  at  low  ambient  densities  could  therefore  be  expected  for  yI  <  2 . 
Taking  7,  =y2  =5/3,  we  find:  (p2  /  p01  )freeze.out  =  1.5  •  KT4 . 

Figure  6(a)  demonstrates  the  time  history  of  a  classical  RM  growth  at 
=  y2  =  5/3  ,  p2  /  p(JI  =  1 0  3 .  Here,  displacement  of  the  contact  interface  8xc  is  shown 


in  units  of  8x0 ,  and  the  normalized  growth  rate  (or,  which  is  the  same,  the  x-velocity 


perturbation  amplitude  at  the  interface)  8vc  -  in  units  of  kut8x0.  Recall  that  typically,  the 


time-dependent  RM  growth  rate  approaches  its  constant  asymptotic  value  quite  fast  [in 
our  present  units  (14),  by  f  =  15  to  20  ],  and  then  stays  near  it,  slightly  oscillating,  see 
Refs.  7,  13,  and  references  therein.  Here,  the  convergence  to  the  asymptotic  value 
apparently  takes  much  longer  —  on  the  time  scale  of  Fig.  6  we  do  not  see  it  yet.  The 
oscillatory  behavior  of  the  growth  rate,  with  two  distinct  frequencies  that  correspond  to 
lateral  sonic  waves  in  two  fluids  separated  by  the  interface,  is  characteristic  of  the  low 
ambient  density  case.  Still,  the  growth  rate  remains  negative,  and  we  observe  an 
approximately  linear  growth  of  8xc  with  some  superimposed  oscillations. 


Figure  6(b)  represents  the  same  for  the  case  p,  /  p01  =  1.5  1 0  4 .  Here,  the  MB 
prescription  for  the  asymptotic  growth  rate  predicts  a  freeze-out,  whereas  the  VMG 
prediction  is  slightly  negative.  Indeed,  the  growth  rate  that  is  initially  negative,  changes 
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sign  at  f  =  1 1 .2  ,  and  the  perturbation  amplitude  reaches  its  peak  value  and  begins  to 


decrease.  Note  that  our  analytic  calculation  has  not  converged  to  the  asymptotic  growth 
rate  (which  must  be  zero  in  the  case  of  freeze-out).  Nevertheless,  comparing  Figs.  6(a) 
and  (b),  we  see  the  growth  rate  shifted  from  negative  to  positive  values  with  decreased 
ratio  p2/p01 ,  so  that  prediction  of  the  freeze-out  seems  to  make  sense.  The  evolution  of 
perturbations  for  the  exact  (that  is,  determined  from  compressible  linear  theory  rather 
than  the  MB  or  VMG  prescriptions)  value  of  (p2  /  p01  )freeze_out  at  T  <  20  should  not 
differ  much  from  that  shown  in  Fig.  6(b).  Note,  that  a  standard-looking  RM  growth  in  the 
negative  direction  that  could  have  been  revealed  by  a  simulation  carried  out  for  the 
conditions  of  Fig.  6(b)  up  to,  say,  f  =  10,  might  be  totally  misleading  —  in  effect,  this  is 
not  a  real  growth  but  rather  just  one  phase  of  decaying  oscillations,  that  should  approach, 
at  very  large  f  ,  some  constant  value  of  5xc . 

III.  NUMERICAL  SIMULATION 

Now  we  compare  our  analytical  results  to  numerical  simulations.  The  comparison 
serves  two  purposes.  First,  we  do  the  comparison  for  a  somewhat  idealized  case  of  a 
constant  intensity  shock  wave  driven  by  external  pressure  into  a  uniform  gas  layer  with 
ripples  at  the  rear  surface.  Here,  the  exact  analytical  results  are  used  to  test  the  accuracy 
of  our  code  in  the  feedout-related  problems.  Second,  we  simulate  the  feedout  produced  in 
a  laser-driven  planar  plastic  target.  In  this  case,  the  roles  are  reversed,  and  the  code  is 
used  to  test  whether  the  predictions  based  on  an  idealized  analytical  model  hold  in  a 
realistic  feedout  situation. 

Our  simulations  were  performed  in  two  dimensions  (2D)  using  the  fast2d 
hydrocode  developed  at  the  Naval  Research  Laboratory14  (more  details  and  further 
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references  are  given  in  Refs.  9,  15).  We  begin  with  a  test  case  of  classical  RM  instability 
developing  at  an  interface  between  two  ideal  gases  with  y,  -  y2  =  5  /  3  and  different 
densities.  The  initial  conditions  correspond  to  the  instant  when  a  shock  front  hits  the 
valleys  of  the  rippled  interface.  The  pre-shock  density  is  p0l  =  1  g/cm3 ,  the  post-shock 
pressure  is  p,=  1  Mbar,  exceeding  the  pre-shock  pressure  by  a  factor  of  1000.  Then  the 
planar  shock  velocity  is  D0  =  1 . 15  •  106  cm/s,  the  post-shock  density  and  speed  of  sound 

are  p{=4  g/cm3  and  ax  =  6.45  •  105  cm/s,  respectively,  as  given  by  Eqs.  (1),  (2)  and 
reproduced  in  the  simulation.  The  wavelength  of  the  interfacial  perturbation  is  A  =  50 
pm,  and  the  initial  amplitude  Sx0  =0.5  pm,  satisfying  the  linearity  condition 
( k8x0  =  0.063  ).  The  time  of  shock-interface  interaction  is  At0  =  2Sx0  /  D0  =0.09  ns. 

Figure  7  compares  the  analytical  and  numerical  time  evolution  of  areal  mass 
perturbation  for  a  conventional  classical  RM  case  of  p2  /  p01  =  10_1  corresponding  to  Fig. 

5(a).  We  see  a  reasonably  good  agreement,  particularly  in  the  times  when  phase  reversal, 
minimum  and  maximum  of  the  mass  perturbation  amplitude  occur  (note  the  early-time 
shift  between  the  two  curves  by  about  A/0 ).  The  difference  between  the  two  curves 
should  be  attributed  to  numerical  diffusion  rather  than  to  a  non-linearity  (even  near  the 
peak  amplitude,  kSxc  does  not  exceed  0.25). 

Figure  8  corresponds  to  the  case  p1  /  p  =  10"3  [Fig.  5(c)].  This  could  be  regarded 

as  a  simulation  of  a  RM  instability  at  large  density  difference,  or,  alternatively,  as  a 
simulation  of  feedout.  From  the  latter  point  of  view,  the  above  small  density  ratio  is 
simply  an  inescapable  numerical  approximation  to  the  actual  (zero)  value.  Dotted  and 
dashed  curves  in  Fig.  8  are  theoretical  predictions  for  p2  -  0  and  p2  /  p  =  10  3 , 
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respectively,  indeed  quite  close  to  each  other.  The  time  of  first  phase  reversal,  t, ,  and  the 


time  interval  between  the  first  and  second  phase  reversals,  At, ,  could  be  compared  to  the 
estimates  (24)  and  (26),  (27): 
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For  7,  =  5/3  we  have:  t,  =  (5  -  V5)/2 ,  £2  =  1/  V2 , 0  =  1.19,  so  t,  =1.7  ns,  At,  =  6.5 
ns,  again,  in  a  good  agreement  with  the  simulation.  The  peak  perturbation  amplitude  is 
reproduced  in  the  simulation  less  accurately  than  in  Fig.  7.  This  appears  to  be  due  to 
greater  influence  of  numerical  diffusion  at  large  density  difference. 

Now  consider  a  60  jtm  thick  solid  plastic  target  with  ripples  at  the  rear  surface 
separating  it  from  a  gas  or  foam  whose  initial  density  p2  is  varied  from  _  to  3/El"6  of 

-5 

the  solid  plastic  density,  p01  =1.07  g/cm  .  We  simulate  irradiation  of  this  target  with  a 
uniform  0.248  nm  KrF  laser  beam,  as  in  Nike  laser  at  the  Naval  Research  Laboratory16. 
The  code  in  this  simulation  includes  inverse  bremsstrahlung  laser  energy  absorption 
(multiple  ray  trace),  Spitzer-H  rm  thermal  conduction,  and  tabulated  equation  of  state 
(EOS)  of  plastic.  Radiation  energy  transport,  which  is  also  included  in  the  code,  has  not 
been  invoked  here,  since  we  are  mostly  interested  in  evolution  of  mass  perturbations  in 
the  bulk  of  the  target,  where  the  radiation  does  not  reach.  We  are  not  attempting  to 
describe  in  any  detail  the  subsequent  RT  growth,  which  could  be  significantly  affected  by 
the  radiation  transport. 

Figure  9  shows  the  results  obtained  for  the  beam  intensity,  which  increases 


linearly  from  1 . 8/E  id  W/cm2to  3/E  I'd  W/cm'  in  6  ns,  and  then  stays  constant  (this 
corresponds  to  an  almost  constant  driving  pressure  Pi  =  0.8  Mbar),  initial  rear  surface 
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ripple  amplitude  8x(j  =  1  pm  and  wavelength  A  =  15  pm.  The  planar  shock  wave 
propagates  at  the  speed  D0  =1.14  •  106  cm/s,  compresses  solid  plastic  by  a  factor  of  2.2, 
and  breaks  out  at  the  valleys  of  the  rippled  rear  surface  at  5. 1  ns.  The  shock-ripple 
interaction  takes  At0  =  2 5x0  /  D0  =0.18  ns,  the  evolution  described  in  Section  IIB  thus 

starting  about  5.3  ns.  The  sound  velocity  in  shock-compressed  plasma  is  al  =1.14  •  106 , 
so  that  the  rarefaction  transit  time  is  2.3  ns,  and  the  acceleration  phase  starts  at  about  7.4 
ns. 

To  check  our  analytical  estimates,  we  have  to  select  some  effective  value  of  y, 
approximating  the  realistic  EOS  of  plastic.  One  can  do  it,  for  instance,  estimating  yt  from 
the  compression  ratio:  according  to  Eq.  (3),  a  2.2-fold  shock  compression  corresponds  to 
Y |  =  8/3  .  Then,  in  the  feedout  limit  p2  — >  0  ,  according  to  Eq.  (37)  (where  7,  =8/3 
corresponds  to  T1  =  1.7  ,  D.  =  0.62 , 6  -  1.77 )  we  obtain:  tx  =  0.36  ns,  A tl  =1.9  ns.  In 
other  words,  the  first  phase  reversal  is  predicted  at  5.6  ns,  and  the  second  —  at  7.5  ns, 
which  is  in  good  agreement  with  the  results  shown  in  Fig.  9.  We  see  that  for 
p,  /  pQ1  <  10"1  the  overshooting  of  the  initial  areal  mass  perturbation  amplitude  in  first 
phase-reversed  peak  is  about  3-fold,  again  in  agreement  with  the  theory  of  Section  IIB. 

At  late  time,  t  — >  °°  the  pressure  p,  would  drive  the  shocked  mass  into  the  half¬ 
space  filled  by  the  gas  of  density  p2  at  some  constant  speed  v2  ,  which  is  found  from  (2)- 
(3),  where  p01  and  v0  should  be  replaced  by  p2  and  v2 ,  respectively.  The  transient 
regime  of  accelerated  motion  lasts  about  (p0I  /  p2)l,2ts  ■  Therefore,  evolution  of  mass 
variation  after  the  rarefaction  breakout  is  determined  by  a  complicated  interplay  of  the 
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RM  instability  at  the  contact  interface,  and  RT  and  RM-like17  instabilities  at  the  ablation 
front.  This  rather  complex  behavior  is  observed  in  Figure  9. 

In  direct-drive  laser  fusion,  a  constant  low-energy  foot  of  the  laser  pulse,  which 
provides  the  desired  value  of  the  target  s  adiabat  during  the  first  shock  compression,  is 
typically  followed  by  a  main  drive  pulse  of  much  higher  intensity.  Most  of  the  target 
acceleration  and  RT  perturbation  growth  occur  during  the  main  pulse.  However,  the  seeds 
for  the  RT  growth,  which  in  our  case  come  from  the  rear-surface  ripple,  are  largely 
generated  earlier.  The  oscillatory  early-time  evolution  of  perturbations  makes  the 
subsequent  RT  growth  phase-sensitive:  a  RT  bubble  can  develop  where  the  target  was 
initially  thinner  or  thicker,  depending  on  the  perturbation  wavelength,  target  thickness 
and  the  laser  beam  parameters. 

This  is  illustrated  by  Fig.  10.  Figure  10(a)  shows  the  laser  pulse  shape,  with  rapid 
intensity  increase  from  the  foot  to  the  main  drive  level  at  about  5.8  ns.  The  background 
gas  density  is  taken  10-3  x  p01 .  In  Fig.  10(b),  the  evolution  of  mass  perturbation 

amplitudes  is  presented  for  several  ripple  wavelengths.  As  expected  [see  (36)],  the  time 
interval  between  the  shock  breakout  and  the  first  phase  reversal  is  proportional  to  the 
perturbation  wavelength  A  .  Shortly  after  this,  the  RT  instability  starts  to  develop,  which 
means  that  the  thickness  of  compressed  target  Ls  becomes  a  significant  scale, 
comparable  to  the  perturbation  wavelength  (in  our  example,  Ls  ~  25  pm).  Indeed,  we 

observe  a  quite  complicated  behavior  rather  than  self-similar.  In  particular,  the 
exponential  RT  growth  occurs  with  initial  perturbation  phase  (albeit  after  two  phase 
reversals)  for  A  =  40  pm,  and  with  inverted  phase  for  A  =  30  pm.  In  other  words,  for 
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A  =  40  pm  the  RT  bubble  develops  where  the  target  was  initially  thinner,  and  for  A  =  30 
pm  -  where  it  was  initially  thicker. 

Figure  1 1  illustrates  the  latter  case  with  the  density  maps  shown  at  t  =  0,  6.9,  8.4, 
and  9.6  ns.  Note  the  first  phase  reversal  at  the  rear  surface  on  Fig.  1 1(b),  similar  to  that 
shown  in  Fig.  1(c).  Then  the  laser  intensity  increase  to  the  main  drive  level  produces  the 
second  shock  wave,  and  by  the  time  the  second  rarefaction  wave  breaks  out  at  the  front 
surface,  the  rear  surface  changes  its  phase  again  [Fig.  1 1(c)].  Although  the  total  mass 
perturbation  at  that  time  has  the  same  phase  as  the  original  rear-surface  perturbation  [Fig. 

1 1(a)],  the  RT  instability  starts  to  develop  at  the  front  surface,  where  the  mass 
perturbation  has  the  opposite  sign,  and  maintains  this  phase  throughout  its  exponential 
growth.  The  3  total  phase  changes  of  mass  perturbation  observed  in  Fig.  10  occur  because 
the  rear  interface  changes  phase  again  and  its  contribution  to  the  total  mass  perturbation 
remains  dominant  before  the  RT  eigenmode  (with  the  opposite  sign)  takes  over.  The  net 
result  is  that  the  bubbles  develop  where  the  target  is  initially  thicker,  in  contrast  with 
what  would  be  expected  in  the  long-wavelength  limit,1’ 3  and  is  actually  seen  in  our  case 
of  A  =  85  pm,  when  there  is  not  enough  time  even  for  a  single  phase  reversal  (or,  to  be 
more  precise,  when  the  phase  reversal  occurring  near  the  rear  surface  is  more  than 
compensated  by  the  RT  growth  starting  near  the  front  surface  at  the  initial  perturbation 
phase.)  Comparing  the  A  =  30  pm  case  to  the  others,  we  note  that  the  two  phase 
reversals  seen  in  Fig.  10  for  A  =  40  pm  occur  because  the  second  shock  reaches  the  rear 
surface  before  the  mass  perturbation  there  changes  sign. 

The  complicated  behavior  of  mass  perturbations  due  to  feedout,  as  described 

2 

above,  has  yet  to  be  studied  in  the  experiments.  So  far,  the  indirect-drive  experiment 


25 


-5 

demonstrated  a  complete  phase  reversal,  whereas  in  the  direct-drive  experiment  the 
onset  of  an  areal-mass  oscillation  at  the  fundamental  harmonic  is  clearly  seen  on  their 
Fig.  4.  To  advance  from  there,  one  needs  to  do  a  careful  design  of  an  experiment  under  a 
number  of  constraints.  1)  In  order  to  observe  more  of  the  oscillations,  it  is  desirable  to 
have  a  longer  driving  laser  pulse  than  those  of  Refs.  2,  3  (about  2.2  ns).  2)  The  laser  drive 
should  be  sufficiently  uniform,  to  prevent  the  target  from  being  torn  apart  by  the  short- 
wavelength  laser-imprinted  perturbations  while  we  try  to  follow  the  perturbation 
development  at  the  fundamental  harmonic  of  the  rear  surface  ripple.  3)  To  increase  the 
frequency  of  oscillations,  it  is  preferable  to  deal  with  short  ripple  wavelength  in  the  range 
15  to  30  pm  (although  by  varying  the  pulse  duration  and  intensity  and  target  thickness, 
one  could  probably  use  longer  wavelengths),  which  is  quite  a  challenge  for  the  imaging 
diagnostics.  In  our  opinion,  the  Nike  laser  facility16  at  the  Naval  Research  Laboratory 
(NRL),  with  its  long  pulses  (4  ns  standard  could  be  extended  to  8  ns)  and  low  beam  non¬ 
uniformity  (0.25%)  is  clearly  the  best  laser  driver  for  such  an  experiment.  The 
monochromatic  x-ray  imaging  diagnostic  technique  developed  and  fielded  on  Nike18 
provides  sufficient  spatial  resolution  to  observe  peak-to-valley  mass  variation  above  2 
pm  of  solid  plastic  at  the  wavelengths  greater  than  1 5  pm.  Design  of  a  Nike  feedout 
experiment  is  now  in  progress. 

IV.  CONCLUSIONS 

We  have  studied  the  feedout  mechanism  that  generates  the  seed  for  the  RT 
instability  starting  from  the  ripples  at  the  rear  surface  of  the  target  through  hydrodynamic 
evolution  of  mass  perturbations  in  a  rippled  reflected  rarefaction  wave.  We  have  shown 
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that  the  feedout  represents  a  quite  special  (stable)  limit  of  a  RM-unstable  flow,  which  is 
only  approached  in  numerical  simulations  if  the  ambient  fluid  density  behind  the  rear 
surface  is  extremely  low.  An  accurate  simulation  of  feedout  requires  much  care  to 
suppress  the  contribution  of  a  numerically  generated  RM  instability  due  to  the  ambient 
gas,  which  is  not  really  there.  The  special  case  of  feedout,  on  the  other  hand,  turns  out  to 
be  analytically  tractable,  so  that  explicit  analytic  formulas  are  now  available  for  the  time 
history  of  the  areal  mass  perturbation  amplitude,  its  early-  and  late-time  asymptotics, 
oscillation  frequency  and  time  intervals  characteristic  of  phase  reversals.  These  formulas 
are  shown  to  agree  with  our  simulation  results.  Since  the  wavelengths  characteristic  of 
the  inner  surface  roughness  in  laser  pellets  are  ~  1 00  pin  or  more,  and  a  typical  thickness 
of  the  shock-compressed  target  is  50  pm  or  less,  the  dimensionless  time  interval 
T,  relevant  for  feedout  in  laser  fusion  conditions  does  not  exceed  ~3  to  5.  Therefore,  most 
of  the  relevant  feedout  processes  occur  before  or  within  the  first  inverted  peak  of  areal 
mass  perturbation.  Amplitude,  timing  and  duration  of  this  peak  are  given  by  our  theory, 
which  agrees  with  our  simulations. 

The  available  data  on  feedout  perturbation  growth  obtained  in  indirect-2  and 
direct-drive3  experiments  agrees  with  the  physical  picture  outlined  in  our  paper.  Much 
more  remains  to  be  done  in  experiment,  however,  in  order  to  study  the  feedout  process  in 
full  detail  and  ensure  validation  of  the  codes  in  the  feedout  regime.  Our  preliminary 
estimates  indicate  that  such  experiments  are  not  beyond  the  reach  of  the  existing  laser 
facilities  and  diagnostic  techniques. 
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Figure  captions 


Fig.  1.  Formation  of  a  rippled  rarefaction  wave  in  a  feedout  situation,  (a)  A  planar  shock 
wave  S  driven  by  a  constant  force  F  approaches  the  valleys  of  the  rear  surface  ripples. 
Mass  and  shock  velocities  v0  and  D0  are  shown  in  laboratory  reference  frame,  (b)  A 
rippled  rarefaction  wave  starts  to  form  near  the  valleys.  In  the  reference  frame  of  shocked 
gas,  its  leading  edge  L  propagates  to  the  left  at  local  speed  of  sound,  ax .  Its  trailing  edge, 
the  expansion  front  E,  propagates  to  the  right  at  constant  velocity  ue .  (c)  A  fully  formed 
rippled  rarefaction  wave  approaches  the  front  surface.  Before  its  breakout,  the  ripple 
amplitudes  at  its  leading  and  trailing  edges,  L  and  E,  are  constant  in  time. 

Fig.  2.  Time  history  of  mass  variation  in  a  weak  rarefaction  wave  at  Mx  =  0.9946  :  exact 
[Eq.  (19),  solid  line]  and  explicit  approximate  [Eq.(20),  dotted  line]  results  of  linear 
theory  normalized  to  pxSx*  ■ 

Fig.  3.  Time  history  of  mass  variation  in  a  strong  rarefaction  wave  at  M \  =  0 .  (a)  Exact 
solution  (22)  and  asymptotic  solution  (25)  (solid  and  dotted  lines,  respectively)  for 
yx  =5/3.  (b)  Exact  solutions  (22)  for  yx  varied  between  7/5  and  9/5.  (c)  Same  for  yx 
varied  between  5/2  and  5:  no  second  phase  reversal  for  yx  >  2.82  . 

Fig.  4.  Time  history  of  mass  variation  in  a  rarefaction  waves  at  Mx  varied  from  0  (thick 
solid  line)  to  0.8  (dashed  line)  with  increment  0.2. 

Fig.  5.  Time  history  of  mass  perturbation  growth  in  a  RM  unstable  configuration:  total 
mass  variation  Sin  and  contributions  to  it  from  perturbations  at  the  contact  interface, 

Smc ,  and  in  the  reflected  rarefaction  wave,  Smr  (solid,  dashed,  and  dotted  lines, 

respectively):  (a)  p7  /  p  =10  1 ,  Mach  number  of  the  reflected  rarefaction  wave 
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Mx  =  0.75  ;  (b)  p2  /  p01  =  10  2 ,  Mx  =  0.53  ;  (c)  p2  /  p01  =  10  3 ,  M,  =  0.36  ;  and  (d) 


Pi  /  Pol  =  10  4 ,  Mx  =  0.24  .  Mass  perturbation  amplitude  is  normalized  to  its  pre¬ 
interaction  value,  (p01  -  p2)Sx 0 . 

Fig.  6.  Displacement  of  the  contact  interface  (  in  units  of  8x0 ,  solid  lines)  and  its  time 
derivative  (the  growth  rate,  shown  in  units  of  knt5x0,  dashed  lines)  for  (a)  p2  /  p01  =  10 
and  (b)  p2  /  p01  =1.5- 10  4 .  Dashed  straight  lines  show  the  asymptotic  growth  rates 
predicted  by  the  MB  and  VMG  prescriptions. 

Fig.  7.  Mass  perturbation  growth  due  to  the  RM  instability  at  density  ratio  p2  /  p01  =0.1. 
Solid  line  —  numerical  simulation,  dotted  line  —  analytical  theory  [same  as  solid  line  in 
Fig.  5(a)]  . 

Fig.  8.  Mass  perturbation  growth  due  to  the  RM  instability  at  density  ratio  p2  /  p01  =  10 
Solid  line  —  numerical  simulation,  dotted  line  —  analytical  theory  [same  as  solid  line  in 
Fig.  5(c)],  dashed  line  —  analytical  theory  for  the  case  of  feedout  [same  as  solid  line  in 
Fig.  3(a)]. 

Fig.  9.  Evolution  of  mis  mass  variation  in  a  laser-driven  60-pm  solid  plastic  target  with 
initial  ripple  amplitude  1  jim,  and  varied  ambient  density  behind  its  rippled  rear  surface 
[here  pback  stands  for  p2  /(1 .07  g/cm  )  ]. 

Fig.  10.  Laser  pulse  shape  with  intensity  increased  from  a  constant  low-intensity  foot  to 
the  main  drive  (a)  and  the  corresponding  evolution  of  rms  mass  variation  in  a  60-pm 
solid  plastic  target  with  initial  ripple  amplitude  1  pm  ,  p,  /  p01  =  1 0  3 ,  and  varied  ripple 
wavelength  (b). 
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Fig.  11.  Density  contours  for  the  case  A  =  30  jim  of  Fig.  10,  taken  at  t  =  0  (a),  6.9  ns 
(b),  8.4  ns  (c),  and  9.6  ns  (d). 
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